function distance = compare_with_standards(I,standard_radii,upper_eyelid,lower_eyelid)
    % Detect the iris
    [center,radius] = circularhough(I,.5,false);
    %eyelid_points = detect_eyelid(I,[center(2) center(1)],radius,false);
    
    [px, py] = circlepoints(radius);
    
    % Scale everything to the radius
    ratio = radius/standard_radii;
    upper_eyelid = upper_eyelid*ratio;
    lower_eyelid = lower_eyelid*ratio;
    
    upper_eyelid(:,1) = upper_eyelid(:,1) + center(2);
    upper_eyelid(:,2) = upper_eyelid(:,2) + center(1);
    lower_eyelid(:,1) = lower_eyelid(:,1) + center(2);
    lower_eyelid(:,2) = lower_eyelid(:,2) + center(1);
    
    % Draw the standards on the iris
    figure; imshow(I); hold on;
    [maxr,~] = size(upper_eyelid);
    
    plot(px+center(2),py+center(1));
    upper_range = upper_eyelid(1,1):.1:upper_eyelid(maxr,1);
    upper_yy = spline(upper_eyelid(:,1),upper_eyelid(:,2),upper_range);

    lower_range = lower_eyelid(1,1):.1:lower_eyelid(maxr,1);
    lower_yy = spline(lower_eyelid(:,1),lower_eyelid(:,2),lower_range);

    plot(upper_range,upper_yy,'r');
    plot(lower_range,lower_yy,'r');
end